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Using a canonical formalism, we determine the equation of state and the phase diagram of eight- 
flavour QCD, as a function of temperature and isospin density. Two mechanisms are at work: 
Bose condensation of pions at high density, and deconfinement at high temperature. We study 
their interplay and find that on our small and coarse lattice the first order deconfinement transition 
appears to end at a critical point at finite density. We investigate the strength of the overlap and of 
the sign problems and discuss implications for the baryonic density case. 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



On the phase diagram of QCD at finite isospin density 



Urs Wenger 



1. Introduction 

QCD at finite isospin density is interesting for several reasons. First of all it is a special case 
of the physically relevant situation where \i u / Hd ^ /As- Secondly, it serves as a platform to assess 
indirectly the limitations of the various numerical approaches to finite baryon density. And thirdly 
it provides a rich and exciting range of physical phenomena: at low temperature and density there 
is a pion gas which transforms into a quark gas at high temperature via a deconfinement phase 
transition, while at low temperature and high density there is Bose condensation of charged pions. 

The transition at finite isospin density is, just like in the case of finite baryonic density, due 
to a large density of a conserved charge, namely the isospin charge Q. However, the system does 
not carry baryon number, but instead the chemical potentials of quarks u and d are set equal in 
magnitude, but opposite in sign, i.e. jl u = —\l c \ = jXj/2. Furthermore, taking the u and d quarks 
masses equal, the system is accessible through lattice simulations because the positivity of the 
theory is guaranteed by 

TirsOysTj =D\ (1.1) 

Using this positivity and QCD inequalities, Ref. [|]] showed that the symmetry breaking related 
to the density driven transition must be accompanied by a condensate (iffiysTi^Y), ie. by %~ ~ 
uy^d, 7t + ~ dy^u states. On the lattice, however, this interesting physical system has largely been 
ignored, with the notable exception of a series of investigations by Kogut and Sinclair [0,BLffl]. 

2. Expectations 

At small isospin densities one can use chiral perturbation theory to describe the physics of the 
system. The chiral Lagrangian is given by 

Se = ^f K Tr[%L%I? - 2/^ReZ] (2. 1) 

where £ G SU(2) is the matrix pion field. The isospin chemical potential /i/T 3 can be included in 
Sfo at leading order without introducing additional low energy constants and breaks SU(2)^ + « — > 
U(1)l+k. The effective potential from the Lagrangian above can be minimised as a function of Hj 
using 

£ = cos a + cos <p + T2 sin <f> ) sin a (2.2) 

where <p is an irrelevant flavour rotation angle corresponding to the residual U(1)l + k symmetry 
and a parametrises the rotation from the standard chiral condensate (ilu + dd) = 2(\j/\j/)o into the 
pion condensate (dy^u) ^ 0. 

From the effective field theory description we can identify two fundamentally different regimes. 
For <m % no pions can be excited from the vacuum and we have £ = 1, i.e. {iiu+dd} = 2(\j/\j/)o. 
This regime corresponds to the standard QCD vacuum at low temperature and zero density. 

For Hi > m K , n + particles can be excited from the vacuum and a Bose condensate of 7T + forms 
where (dy^u) ^ 0. In that regime the remaining U(1)l + r symmetry is broken spontaneously and 
hence the transition belongs to the universality class of the 3d XF-model if it is second order. 

Finally, for Hj > m p the chiral expansion breaks down and the effective field theory descrip- 
tion eq. ( |2.1| ) is no longer applicable. Nevertheless, on general ground one can argue [[!]] that the 
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Figure 1: The free energy difference (F(Q) — F(Q — l))/T versus the number of M-quarks Q. Left: for a 
gas of free (blue) and weakly interacting (black) pions and for a gas of free massless fermions (red). Right: 
simulation results on 8 3 x 4 lattices for successive temperatures, increasing from ~ -f to ~ 1 . IT C from top 
to bottom. 

condensate remains {dy^u) ^ fj and one might expect a BEC-BCS crossover at large jj,] when the 
pions "dissociate" (the QCD running coupling becomes weak). 

A complete description of the various phases is provided by the equation of state (EoS) of the 
theory. Here we do not attempt to calculate the full EoS but rather the relation between the isospin 
density and the chemical potential 1 , 

P/ = ^= P/ ( M/ ). (2.3) 
The canonical free energy F(Q) = — \nZc{Q) and its derivative 

F(Q)-F(Q-l) V =^ ^~=Vi (2-4) 

dpi 

are numerically well accessible from a canonical formulation of lattice QCD. 

The behaviour of the system in the various expected phases can be anticipated using simple 
semiquantitative considerations. For example, in the confined regime at low temperature and low 
density, we have a gas of (to a first approximation) free pions. For bosons of mass m = mT and a 
chemical potential /J. = fiT, the density reads 

Pfree(A,-) =^j*PP (j^f^i ~ ^T^~[ ) ^ 



where CO = \f p 2 + rti 2 . The behaviour at zero and low temperature is sketched in Fig. [I] where we 



plot the free energy difference, eq. ( ]2.4[ ), versus the number of w-quarks Q. At /I/ = m n the system 
Bose condenses and as a consequence the free energy difference remains constant beyond that point 
(denoted by the circle in Fig. [[]). For an interacting gas of pions the 'EoS' can be described by the 



effective field theory eq. (EJh and the density at zero temperature is given by 



P/ =/^/(i-(^r). (2.6) 



In a slight misuse of notation we will nevertheless refer to eq. (2.3) as EoS in the following. 
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The interaction leads to an increase of the free energy difference at a given density and temperature 
as sketched in Fig. [j} 

In the high temperature deconfined phase and for quark mass m q <C £t/ the pressure can be 
approximated by that of a massless, free Fermi gas, leading to 

§ = £/ + ^£/ 3 - (2-7) 
The corresponding situation is sketched in Fig. 0. 



3. Numerical results 



For computational convenience we consider in our numerical simulations two staggered fermion 
fields corresponding to Nf = 8 QCD in the continuum, i.e. four u- and four d— quark species, all 
degenerate in mass. We keep the lattice volume fixed to 8 3 x 4 as well as the bare quark mass 
am = 0.14 while varying the bare coupling j3 and the isospin chemical potential a\i\. Note that the 
bare quark mass is tuned in such a way that the deconfinement transition at T c and zero density is 
first order. Varying j3 changes the lattice spacing a and thereby the temperature T = i and we 
cover a range of temperatures between \T C <T < T c . We emphasise that for the range of couplings 
we consider, the pion mass am % , measured at zero temperature, changes very little in units of the 
lattice spacing, so that the ratio nix/T stays almost constant. Note that, contrary to Refs. [§, ||, |J, 
we do not introduce a U (l)-breaking source term in our action. Thus, we can observe the spon- 
taneous symmetry breaking characteristic of Bose condensation without performing any delicate 
extrapolation. Finally, our results are obtained by combining 68 ensembles at six values of jii up 
to jJ-i/T < 5 with Ferrenberg-Swendsen reweighting. 

The free energy is obtained from the canonical partition function, F(Q) = —\/T\ogZc{Q). 
Zc(Q) is estimated using standard grand canonical Monte Carlo simulations 



Z C (Q) _/ |det e | 2 



(3.1) 



Z GC {]Xi) \|detDGu/)| 2 /z GC 

where detg is obtained by decomposing exactly the part of the measure which depends on \ii 
i.e., the fermion determinant, 

+3V 

detD(/x 7 )= £ det e exp(g/i//r) . (3.2) 

Q=-3V 



Our results for the free energy difference eq. < J2-4| > (in units of T) versus the number of w-quarks 



Q for successive temperatures together with the theoretical expectations are shown in Fig. [I] right. 
At high temperature (lowest curve) our data is very well described by an ansatz motivated by 



the free Fermi gas description, eq. ( |2.7[ ), where the coefficients for the linear and cubic term are 
multiplied by factors that take into account effects of interactions. The numerical values we obtain 
are reasonably close to one and consistent with those obtained in [^]. 

Decreasing the temperature successively, we obtain the second-, third-, ... to-lowest curves. 
The S-shape of the curves for T <T C indicates that P/(/i/) is multivalued, and is thus characteristic 
of a first order transition - in the thermodynamic limit the density would jump at the critical isospin 
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Figure 2: Phase diagram in the (fa/T,T /T c ) -plane from grand-canonical Monte-Carlo simulations (GCMC) 
at finite isospin chemical potential. The red circle indicates the critical point as inferred from the change 
of the behaviour of the Binder cumulants (shown in the inset for the quark condensate) at fa/T ~ 2.5. The 
green circle indicates the transition into the pion condensed phase (BEC) as obtained from the canonical 
analysis. 



chemical potential. Here, the critical fa can be determined using the Maxwell construction and 
it increases as the temperature is decreased. Furthermore, the S-shape becomes smoother and 
finally disappears at around T /T c ~ 0.8 where the jump in the density vanishes and the first order 
transition is replaced by a crossover. In order to ascertain the presence of this critical point, we also 
monitor the distributions of the plaquette, the Polyakov loop, the quark condensate and density by 
measuring the corresponding Binder cumulants in grand-canonical simulations at fixed fa. These 
cumulants remain essentially constant in the region < fa/T < 2 beyond which they start to grow 
linearly, cf. inset in Fig. ||| Attributing this behaviour to the critical point of the 3d Ising model 
universality class one can locate the point at fa/T ~ 2.5. A finite-size scaling analysis and a 
continuum extrapolation are of course required to locate this critical point reliably, in particular, 
to disentangle it from the 2nd order BEC transition. Our findings are similar to those of Ref. [[|], 
which for Nf = 3 also find that the Binder cumulant grows as an isospin chemical potential is on. 

Lowering the temperature further, we observe a transition from the pion gas into the Bose 
condensed phase, marked by a circle in Fig. [I] (right). It turns out that the system can very well be 
described by the ansatz in eq. (2^) where the value of m % is fixed from our measurement at zero 
temperature and f n is treated as a free parameter. 

The full phase diagram in the (fa/T,T /T c ) -plane is presented in Fig. || where the black circles 
denote our results from the grand-canonical Monte-Carlo simulations (GCMC) and the red circle 
marks the critical point. For fa < m n the pseudo-critical line between the hadronic gas and the 
plasma can be well described by a quartic fit (dashed line). The green point marks the entry into 
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Figure 3: Left: Pion susceptibility Xx+ f° r various values of Right: Data collapse obtained from shifting 
j8 to P c and rescaling both the x— and axis. Also shown is the universal function from simulations of the 
3d XT-model. 

the pion condensation phase obtained from Fig. [I] (right). As expected, for T > this point is shifted 
towards larger values of Hi as compared to the zero temperature situation where the transition would 
occur exactly at Hi/T = m n /T (indicated by the small arrow on the x-axis). 

The transition from the plasma to the pion condensation phase can be determined by moni- 
toring the pion susceptibility % n + = L, T (7T + (0)7r + (x)). Since this transition is characterised by the 
breaking of a U(l) symmetry, it is expected to belong to the universality class of the 3d XF-model. 
We can expose the universal behaviour of the transition by comparing the data at different values 
of /I/ with the universal scaling curve obtained from simulations of the 3d AT-model. Fig. || shows 
that the agreement is excellent. 

Finally, our calculation also provides a crosscheck of the reweighting method from zero to 
finite (isospin or baryonic) chemical potential. It turns out that for isospin chemical potential 
reweighting from pLi = is possible, but fails to describe the phase diagram reliably. The overlap 
of the /i/ = ensembles with the ones relevant at finite density is not sufficient to determine e.g. the 
order of the phase transition at finite jU/, and the extrapolation to finite pLi noticeably underestimates 
the density of the system. This is understandable from the fact that one tries to extract finite density 
information from the tail of a density distribution which is generated at /I/ = 0, i.e. centered at zero. 

For an extrapolation to finite baryonic density one has the additional problem that the average 
sign of the fermionic determinant tends to zero very quickly with jj,, as illustrated in Fig. hence 
narrowing the range of applicability of the reweighting procedure even further. 

4. Conclusions 

We determined the EoS and the phase diagram of Nf = 4 + 4 QCD at finite isospin density 
and finite temperature. The two mechanisms at work are clearly exposed: we observe Bose con- 
densation of pions at high density and deconfinement at high temperature. For our quark mass, 
the transition from the hadronic gas to the quark-gluon plasma is first order at zero density. As the 
density is increased, the transition appears to turn into a crossover at pii/T ~ 2.5, in qualitative sim- 
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Figure 4: Left: same as Fig. 1 right, but reweighted from zero density simulations only. Right: isocontour 
lines of the average sign (detg /|detg| 2 }. 

ilarity with the results of Ref. []3p. We caution, however, that to confirm these findings a finite-size 
scaling analysis is still needed. 

The transition from the hadronic gas to the BEC phase is determined from a canonical analysis 
of the free energy density, occurs at low temperature T/T c < 0.5 at /»/ slightly larger than m % and 
is consistent with being second order. The transition from the BEC phase to the plasma phase is 
also second order and its universal behaviour is confirmed to be that of the 3d XF-model. 

Finally, by comparing our (interpolated) results for the EoS P/(/i/) with the ones obtained 
by reweighting from zero to finite jii, we find that the latter technique systematically underesti- 
mates the densities as the chemical potential grows, and that to determine its range of reliability a 
comparison of different approaches is needed. 
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